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Abstract 

We report a statistical analysis of the Island ECN (NASDAQ) order book. We de- 
termine the static and dynamic properties of this system, and then analyze them 
from a physicist's viewpoint using an equivalent particle system obtained by treat- 
ing orders as massive particles and price as position. We identify the fundamental 
dynamical processes, test existing particles models of such markets against our find- 
ings, and introduce a new model of limit order markets. 
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1 Introduction 



Physicists can contribute much to the study of interacting humans, because of 
their familiarity with cooperative phenomena. Collective economic behavior is 
an area in which physicists have recently become acutely interested, probably 
because of the availability of a huge amount of free electronic data related 
to economy. In particular, physicists' attention [1] has focused on financial 
markets [2-4]. 

Most of the time, the only information one can obtain about a particular 
stock is the price and volume of the last transaction, and the current price 
and number of shares of the best buying offer (also called the highest bid), 
and, of course, the price and number of shares of the best selling offer (lowest 
ask). 

However, in electronic limit order book marketsQ, more data are available, 
either publicly or by request: one knows not only the highest bid (buy) and 
the lowest ask (sell) orders, but also several lower bid and higher ask orders. 
This obviously gives more information than only the last price and volume. 



^ Like for instance the NASDAQ, the Bourse de Paris, the Swiss Market (SMI), 
the Australian Stock Exchange (ASX) , . . . 
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but it is several fold more difficult to analyze: in addition to the last price and 
volume time series, one has to deal with two time-evolving unidimensional 
order distributions. 

There is growing literature on the analysis of limit order book markets, partic- 
ularly in Economics [5]. In Econophysics however, this area is almost a terra 
incognita, except for a very recent work [8], where the five best prices on both 
bid and ask sides of the global NASDAQ book are analyzed. 

Here, we analyze data obtained from a sub-part of the NASDAQ, namely the 
Island Electronic Communication Network (ECN) (www.island.com), which 
displays in real time the 15 best orders of both typesj^]. While previous studies 
have mainly addressed static properties of such markets, we provide a dynamic 
analysis that is able to identify some fundamental properties of such markets. 
In this way, by analyzing the available data from an appropriate point of view, 
we are able to arrive at the necessary ingredients for a minimal market model. 



2 Definition of the market 

There is one order book per stock in a given ECN (but there are several ECNs 
that are part of the NASDAQ); it reflects the current state of investors public 
wishes by listing what quantity of shares they are willing to buy or sell and 
at which price. An order k is characterized by three intrisic quantities: 

• its type, bid (buy) or ask (sell); for each order k this can be represented by 
a binary variable G {B, A} = {+1, —1}; 

• its price pk, which is discretized; 

• its size mfc, in shares. 

In addition, the owner of an order can specify a timeout, time interval after 
which the order is automatically removed if the order has not been ffiled. The 
life of this order {akiPki^k) begins at time t\ when it is placed at price pk- 
Possible prices are discretized, and the elementary step is called the tick — 
in the Island ECN, a tick is equal to $ 1/256. It appears on the order book 
until time t% when it is removecQ, or when a new order {—ak,Pk,fn') of the 
opposite type and of at least the same size (m' > mk) is placed at the same 
price, leading to a transaction!^; order k's lifetime is then = t% — t\. Note 
that about 9 ± 1% of all orders happen to be partly filled at some point, that 
is, since t\, an opposite order of a smaller size has been placed at price p^. 

^ Note that two orders of the same kind may have the same price. 

^ This may be due to a cancellation by its owner, or by its automatic removal, once 

its time-out is reached. 

^ This is called a price-taking behavior. 
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Therefore, the size of an order is not a conserved quantity. There are of course 
transaction costs: Island ECN charges $ 0.0025 per share for a price taking 
orders, and $ 0.0015 per share for the other ordersQ, but removing an order 
is free. 



3 The data 



We collected high frequency!^ data from the Island ECN order book of four 
stocks, namely Cisco (CSCO), Dell (DELL), Microsoft (MSFT) and Worldcom 
(WCOM). Before doing any analysis, we want to warn the reader that almost 
all results of our study may be affected by the fact that only the 15th highest 
bids and lowest asks orders are displayed at a given time. We checked that 
most of the activity takes place within the range of the 15th best orders, and 
we managed to keep track of consistently more than 75% of all orders, but 
the reconstructive method we designed distorts the reality to an extend that 
is hard to estimate. Finally, note that we may not pick up an order whose size 
is greater than 10000 shares since it does not need to be displayed, even if it 
exists in the market. 



4 Static properties 



An order is anonymous, i.e. there is no way to identify its owner, but each 
order has a serial number, allowing us to study the fate of every order. 



4-.1 Order size 



It is well known that orders have a tendency to cluster both in size and position 
[5]: they tend to have a size of multiple of 10, 100 or 1000, and to be placed 
at round prices, or at halves. Fig 1 reports P{m), the probability distribution 
of their size for various days of CSCO, and clearly shows this clustering. We 
found also that the density of selected sizes does not depend on the stock nor 
on the day. 
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Fig. 1. Order size probability distribution P(m) of the bid size of DELL (10.04.2001). 
The clustering is obvious. 




1 10 100 1000 10000 1 10 100 1000 10000 



Fig. 2. Lifetime histogram for bid (x) and ask orders (circles) of all orders of all 
stocks and all days. The left graph plots -ffev, the histogram for cancelled orders 
(straightline: power-law with a —2.1 exponent). The right graph plots i/aii) tbe 
histogram for fully filled orders (straightline: power-law with a —1.5 exponent). 

4-2 Order lifetime 



Fig 2 reports the histograms of the lifetime of all orders in our data. The 
lifetime distribution for cancelled and fully filled orders seem to have an alge- 
braical decay with exponents that coulcQ be — 2.1 ± 0.1 and — 1.5 ± 0.1 The 
lifetime distribution of cancelled orders has visible peaks at r = 90, 180, 300 

^ thus rewarding patient investors. 

^ Every one or two seconds on average. 

^ Our data actually only provide only an estimation of the lifetime of each order, 
hence the value of these exponents are only indicative. 
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Fig. 3. Histogram of aggregated orders for CSCO on 14.12.2000, showing that the 
orders are separated typically by $ 1/16 although the tick is $ 1/256. 

and 600 secondsQ it is clear that they correspond to typical timeouts. It is 
very likely that another typical timeout is 60 seconds, although this peak does 
not clearly appear on Fig. 2. 

4-3 Distribution shapes 

The shape of the orders distribution is a dynamical quantity in essence and can 
change tremendously in a few seconds. Plotting time-averaged distributions 
reveals some general features. For all time t, we define At{p), and Bt{p), which 
are the number of shares to be sold, and bought, at price p. Fig 3 plots the 
functions a{Ax) = {At^Axas]^)) and b{p) = {Bt{Axud)) , where Ax is the 
price relative to the relevant extremum of the considered distribution, that is, 
Axaskit) = pit) — am(t) > where am is the price of the minimum ask, and 
a;bid(^) = p(t)—bM(t) < 0, where Bm is the price of the maximum bid. Note that 
the so-called bid-ask spread is then simply g(t) = amif) — bMif) The clustering 
in the position of the orders clearly appears on fig 3: although the tick is 
$ 1/256, there are clear peaks separated by 16 ticks — that makes $ 1/16, 
which is the official tick of the NASDAQ. Consequently, the effective tick is 
$ 1/16, and one has to coarse grain the raw distributions a{Ax) — > a' [Ax) 
and b{Ax) b'{Ax) by a factor 16; this is done on fig. 4. 

It appears that the distributions are convex and peaked at Ax = 0. This could 
explain the origin the overdiffusive behavior of prices which has been observed 
in various financial markets [3,4]. Indeed, suppose that the order distributions 
are fixed during a time window and that the price follows a random walk, 

^ The nature of our data broadens these peaks. 



5 



CSCO 14.12.2001 



MSFT 26.10.2000 



1500 




Fig. 4. Histogram of aggregated orders of CSCO on 14.12.2000 (left graph), and 
MSFT on 26.10.2000 (right graph). The left distribution is the buy side and the 
right one the sell side. 

then the convexity of the orders distributions causes an over-diffusive price 
behavior. 

We could not find a general mathematical shape for distributions a'{Ax) and 
b'{Ax), maybe due to the fact that the data we have is not reliable for the tails 
of the distributions; we could note however that the buy and sell distributions 
are usually not symmetric, and this shows up through the virtual impact 
function discussed in the next subsection. 



4-4 Market Impact 



Market impact is an important quantity: it is defined as the price change 
Ap caused by a large order[^ of size M, hence, it is a functional relationship 
between the size and the price change. There are actually two types of impact. 
The first one is the virtual impact, which is related to the integral of the order 
distributions: Ap{M) is defined through M = \ r{x) dx\, where r{x) = 
a{x) or b{x). Since apparently the average shape of the distributions have no 
general characteristic form, one can expect that the virtual market impact 
function differs from stock to stock and from day to day. This is the case, as 
shown in fig. 5. In the literature, it has been argued [14] and suggested from 
data [8] that this impact should have a power-law form, with an 0.5 exponent 
[14] , or with an exponent near 2 [8] . From our data, which are not reliable for 
the distribution tails, we could infer that the virtual impact functional form 
may not be strictly speaking a power-law, but can be approximated by such a 
form, with an exponent ranging from 1 to 3, depending on the day and stock. 

^ Usually, large orders are split into several smaller orders in order to reduce the 
market impact. 
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Fig. 5. Virtual impact of four days of CSCO. Circles correspond to the ask side, and 
diamonds to the bid size 

On the other hand, the real impact is defined as the real price change for a 
real large order, or a fast succession of almost simultaneous orders; the two 
types of impact can significantly differ. 



5 Dynamical properties 

Dynamical analysisp"^ is made easier by a straightforward analogy between 
orders and particles: an order can be thought as a particle whose mass is the 
order's size, and whose price is a spatial position. Since the possible prices 
are discretized, it makes sense to consider particles on a lattice whose mesh 
size is equal to the real tick. Placing an order is equivalent to depositing a 
particle, withdrawing an order amounts to evaporating it. A transaction takes 
places whenever two particles of opposite type are deposited at the same price, 
most often sequentially: this can be thought as a annihilation; note that an 
annihilation is caused by a deposition, hence there are two types of deposition 
depending on where an order is placed. 

This analogy is already used for the modelling of limit orders markets in [9,10], 
where two types of particles-orders diffuse on a lattice and annihilate whenever 
they meet, and in [12], where orders do not diffuse, but are randomly deposited 
and annihilated. Hereafter, we consider sell (ask) orders as type A particles, 
and buy (bid) orders as type B particles. The terminology used hereafter is 
that of physics; in particular, we will use the word gap instead of spread. 



The analysis presented here can be compared with that in ref [6]. 
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Fig. 6. On the left, from top to bottom: bid rates of deposition 6, evaporation rj and 
annihilation a in events /second, time averaged over a window of 1500 seconds of a 
day of DELL. Right: ratios 6/r] and 6/a, and 6 — r] — a 



5. 1 Temporal properties; rates 



We consider three types of event, namely deposition, evaporation and annihi- 
lation]^, and denote their rate by (5, r] and a respectively. Whenever needed, 
subscript B labels rates associated with S-type particles, and A with A-type 
particles. A rate is defined as the average number of events per second. 

The method we used to plot fig 6 was to count how many events occurred 
during a given time window, typically Dt = 30 seconds, then of course to 
divide this number by Dt, and to run a moving average in order to smooth 
the graphs. Figure 6 reports the B rates of one day of DELL. They fiuctuate 
wildly during the day, but not in an independent way: the autocorrelation 
and crosscorrelation functions of S, rj and a are slowly decreasing, apparently 
algebraically, as shown by fig. 8 for 6, rj and a; the gap was also found to share 
this property. It is tempting to consider these facts as the origin of the well- 
known algebraical decay of the autocorrelation of the price volatility. Note 
that {6{t)r]{t - At)) > {5{t)r]{t + At)) {At > 0), which is due in part to the 
redeposition of evaporated orders (see also section 5.4). The distribution of the 
number n{t) of events during t and At has exponential, hence Poissonian tails 
(fig. 7), but this distribution is computed over a whole day, and is therefore 
a superposition of a time-varying distribution. The fact that a Poissonian 
distribution is obtained is merely a consequence of this superposition and the 
associated loss of information. 

We emphasize that the evaporation rate t] is high, which is consistent with the 



An annihilation is always caused by a deposition of an opposite particle, but we 
do not include this type of event into the deposition rates. 
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Fig. 7. Cumulative probability of event rates 5 (circles), r] (squares) and a (dia- 
monds) of bid (filled symbols) and ask (plain symbols) sides (DELL 09.04.2001) 
{Dt = 30 s). 




Fig. 8. Bottom graph: autocorrelation of the bid rates 6 (top), r] (middle) and 
a (bottom) measured with At = 1 s the continuous line has a -0.17 exponent. 
Upper graph: cross-correlations {5{t)r]{t + At)) (top) and {6{t)a{t + At)) (DELL 
09.04.2001) 

short typical order lifetimes and the fact that removing an order from Island 
ECN is free of chargeP^. We found that typically 80% of all orders evaporate, 
and 20% anihilate, in agreement with results of ref [6], where similar evapo- 
ration rates were also measured in the OM Stockholm and London Securities 
and Derivative Exchange . This implies that the absence of evaporation in 
model of ref. [12] is an important shortcoming. 



Hence the typical evaporation rate may be smaller in ECNs which charge for 
cancelling orders. 
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Fig. 9. Frequency of an event versus its relative position (DELL 09.04.2001) (left 
column: bid side; right column: ask side) 

There are some robust properties of the dynamics of the books: the right 
column of figs. 6 shows that the ratio rj/S is quite constant during the usual 
trading period. On the other hand, the ratio a/ Sis obviously related to the real 
trading activity and shows periods where it remains constant, but can greatly 
vary during a day, suggestting that the composition of the market's population 
has changed. The way the rate a is defined implies that an annihilation does 
not implies a loss of a particle (orders can be partly filled), hence r) + a 
is consistently greater than 5 and the particle conservation law should read 
1]/ + a' = 6, where a' /a = est . On the other hand, the mass rates d^, 1]'" 
and a^, which count how many shares deposit, evaporate and annihilate per 
second, often obey a mass conservation relationship 6" = rj"" + a" during a 
trading day; we checked that these rates also have algebraically decreasing 
auto- and cross-correlation functions. 



5.2 Spatial properties 



It is known that the activity in order driven markets occur near the relevant 
extrema of the distributions [5,7]. The shape of the orders distribution (see 
section 4.3) is somewhat a confirmation of this statement, but the rate of 
activity is an even better one. Fig. 9 reports the frequency of the three fun- 
damental events versus their relative position Ax. It is clear from that figure 
that all events usually occur close to Ax = 0. By comparing the probability 
of an event at position Ax (fig. 9) with the average number of shares at that 
position (fig. 4), which is in first approximation proportional to the number 
of orders, one can see that the probability of an event depends on x. 
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Fig. 10. Upper graph: spatial width of deposition versus gap averaged over all stocks 
and all days (bid: . Lower graph: normalized integral W{g)/g of w{g) (bids: thick 
curve, asks: thin curve), and a reference line y = 0.18 x + C (dotted line). 

5.3 Spatio-temporal properties 

This type of property is much harder to analyse. Nonetheless, it appears that 
the deposition process depends on the gap g. Indeed, defining the spatial 
width of deposition w = ^ (Ax^gp), one obtains fig 10 which reports wa^q) 
and WB^g)- We find that w is constant up to (7 ~ $ l/16,Qf] i.e. as long as 
the two order distributions are adjacent; when g > $ 1/8, i.e. when there is 
an effective tick between the two distributions, we found a linear relationship 
between w and g for small values of g, as shown by W{g)/g = Jq w{g')dg' /g, 
the normalized integration of w{g), which suppresses the noise. 

We found a linear relationship for both orders that are deposited onto the gap, 
and the ones that are placed in the bulk of the distributions. It is easy to see 
how gap orders lead to this relationship: suppose that some trader is willing 
to pay less than the current best offer, but still in a hurry, then she will place 
an order anywhere in the gap, e.g. at the middle point, causing the observed 
effect. On the other hand, it could be possible that the linear dependence 
measured for bulk deposition were due to the small price window defined by 
the 15 best orders on each side; however, we found the this window is typically 
of order $ 0.8, whereas fig 9 shows that most of the activity occurs in a small 
partof the price window. This is an evidence that allows us to conclude that 
this linear dependence is genuine. 



this is the real basic length scale of this system (see above). 
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Fig. 11. Histogram H{D) of individual order diffusion coefficients. The vast majority 
of orders diffuse less tfian one effective tick during their life. 

5.4 Diffusion 

We could also test if orders diffuse on the lattice, i.e. if an order jumps from 
a price to another before meeting an order of opposite type and annihilating, 
or simply evaporates. Orders diffusion is the central assumption of models 
found in ref [9,10]. The method we used is to keep track of orders with un- 
usual sizes, for instance 832 shares. When such an order first evaporates, we 
begin to record its trajectory, until it is annihilated. We found that about 2/3 
of such orders are never deposited again, i.e. they simply vanish; if they are 
deposited again, their average number of redeposition is about 2; amongst the 
redeposited orders, about half of them are redeposited at the same price a few 
seconds later, which shows up in the asymmetry of the cross-correlation func- 
tion {6{t)ri{t — At)) (fig 8). The diffusion coefficient D of the diffusive orders 
turns out to be typically D < 10^^[$^/s] (see Fig 11), which corresponds to 
an average price deviation of one effective tick during 1000s [^. This suggests 
that the fundamental process of the orders dynamics is not diffusion; however 
such an approach is useful, particularly concerning analytical results [9,11,10]. 



6 A new breed of models for limit orders markets 

We believe that the above facts have strong consequences for the modelling of 
limit order markets. They imply that one has to reconsider the two existing 

Our method actually sometimes overestimates the diffusion coefficient, since it 
may happen that a odd-sized order that evaporates belongs to another person when 
it is deposited some time later. 
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families of non-strategic models of limit orders, those of Bak, Paczuski and 
Shubik (BPS) [9,10], the one of Maslov [12]. Both are very appealing thanks 
to their extreme simplicity, which makes them suitable for analytical anal- 
ysi^33 they are able to mimick some typical market behavior, such as 
volatility clustering and large price variation probability density function, but 
both include fundamental hypotheses that are in part incompatible with our 
findings. In addition, at least under their simplest form, they have a price be- 
havior that differs from what is observed in real markets. Indeed, the returns 
^Atit) between time t — At and t, usually defined as logp(t) — log p{t — At), 
where p{t) is the price at time t, are known to be overdiffusive for short times, 
i.e. (rAf)t oc t-'^^^*) with H{At) > 1/2 for small At, and to tend to those of a 
random walk for large times, i.e. H{At) ~ 1/2 for At oc a year [3]. The basic 
BPS model and Maslov's model are characterized by H{At) = 1/4 for all At, 
hence the prices are underdiffusive and do not tend to a random walk in the 
long time limit. However, as soon as the orders diffusion is biased towards the 
gap in BPS model, H = 1/2 is recovered for large times At. Obtaining an 
over- diffusive behavior seems to require an (strategic) imitating behavior that 
leads to some kind of herding [9] . 

Our aim is to build a model compatible with our observations. We retain 
two processes, namely deposition and evaporation: an annihilation being first 
caused by the deposition of a new order on top of a previously existing op- 
posite order, it can be included into the deposition process. In addition, we 
consider that all positions are relative to the relevant extremum of the order 
distributions, i.e. am(t) and 6^/ (t)[^. For the time being, we restrict the mass 
of orders to be all the same, and consider all events to be independent; we 
also suppose that only one order can be deposited during a timestep, that the 
rates are constant, and shall discuss the consequence of each hypothesis in the 
following. 

6. 1 Definition of the model 

We consider an infinite one dimensional lattice and two types of particles, A 
and B corresponding to the types of order. These can be considered as partic- 
ular states of X— particles (X G {A, B}); two types of events can happenp^: 

• Deposition: with probability 5, a X— particle is deposited at site x drawn at 
random from the pdf Px{x,t). Whenever the particle X is deposited into 

See ref [11] for an analytical approach of BPS models, and [13] for a mean-field 
solution of Maslov's model. 

This contrasts with Maslov's model, where events occur relative to the last paid 
price, but may not be crucial. 

All quantities below can of course depend on the type of particle. 



13 



12000 




11000 I ' ' ' ' ' 1 

10000 20000 30000 

time 



Fig. 12. Two time scries of the market model with parameters a = 1/2, S = 1/100, 
K = 2, C = 3 (upper graph) and K = 1 (lower graph), showing the dramatic 
influence of K on the system's properties. 

the bulk of the —X distribution, one considers that the annihilation takes 
place at the relevant extremum of the —X distribution 
• Evaporation: each X— particle on the lattice has a probability rj to evapo- 
rate. 

In order to complete the definition of the model, the pdf Px{x,t) has to be 
specified. For the sake of simphcity, we consider a Gaussian function.. For 

the time being, we assume that it is centered on the relevant extremum of 
the particle distribution. Its variance ax{t) can cither be constant, or linearly 
depend on the gap: ax{t) = Kg{t) +C; if C = 0, fixes the annihilation rate 
a, and if C 7^ 0, a increases as the gap decreases. The fact that every particle 
has the same probability of evaporation implies that the number of particles 
has a well defined average and may be seen as equivalent to a time-out of 
order 1/r]. This model has therefore four parameters, 5, r], K and C. 



6.2 Results 



We observed power tails of the returns and, depending on parameters, volatil- 
ity clustering. The parameters have strong effects on the dynamics of the 
model (see figure 12). As soon as evaporation takes place, the returns have a 
random walk behavior at large times (fig. 13). This means that evaporation is 
probably enough to remedy the underdiffusive behavior of prices even at large 
times of Maslov's model. 

One of the main assumptions of the model we propose is that at most one or- 
der is deposited during a timestep. This is in clear contradiction with markets 
data. Let us extend our model in order to include simultaneous depositions. 
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Fig. 13. Results from the market model: typical (r^t) versus At, showing that the 
Hurst exponent H is 1/4 for small times and tends to 1/2 for large times (circles: at 
most one particle is deposited during t and t + 1 S = 0.5, rj = 0.01, K = 1.5, C = 10; 
x: simultaneous deposition of an average of 5 particles of both types per unit of time 
(same parameters); ★: simultaneous deposition 6 = 2, r] = 0.04, K = 1.5, C = 10. 

If at each time step, the number nx{t) of deposited orders of type X is drawn 
from an exponential distribution with average S, the region where the price is 
under-diffusive becomes gradually less pronounced when 6 is increased, even 
while keeping 6/7] constant. For the time being, the present model has no 
overdiffusive behavior, which seems to be a property shared by all simple mar- 
ket particles models with no strategic behavior. The effect of the non-trivial 
time properties of the real rates is still unclear, but being the consequence of 
agents behavior, it may be an indirect origin of the overdiffusive behavior of 
financial markets. In addition, we note that in our model as well as in Maslov's 
model, the presence of volatility clustering is due to the algebraical decay of 
the gap autocorrelation function. Finally, the effect of heterogeneity in order 
size is still to be determined. 



7 Conclusions 

We have determined several features of Island ECN's order books. Some of our 
findings may be specific to this ECN. The static features we studied include 
the size and life-time probability distributions, the average shape of the orders 
distribution, from which we derived the several market impact functions. More 
importantly, we have been able to determine fundamental mechanisms of such 
markets, which has important consequences for the modelling: this lead us to 
consider alternative over-simple particles models, which offer a closer descrip- 
tion of real markets. We have observed that the under-diffusive behavior of 
simple particles models is less pronounced if the deposition of more than one 
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particle per time step is allowed to deposit. Finally, a non-strategic behavior 
that leads to overdiffusive price behaviors in particles models has still to be 
found, and may require non-constant rates with non-trivial interdependence. 

We thank Y.-C. Zhang, S. Maslov, M. Marsili, and Th. Bochud for useful 
discussions. This work has been supported in part by the Swiss National Funds 
of Scientific Research. 
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